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Abstract. We study analytically and numerically the problem of a nonlinear mechanical oscillator with 
additive noise in the absence of damping. We show that the amplitude, the velocity and the energy of 
the oscillator grow algebraically with time. For Gaussian white noise, an analytical expression for the 
probability distribution function of the energy is obtained in the long-time limit. In the case of colored, 
Ornstein-Uhlenbeck noise, a self-consistent calculation leads to (different) anomalous diffusion exponents. 
Dimensional analysis yields the qualitative behavior of the prefactors (generalized diffusion constants) as 
a function of the correlation time. 

PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.10.Gg 
Stochastic analysis methods (Fokker- Planck, Langevin, etc.) - 05.45.-a Nonlinear dynamics and nonlinear 
dynamical systems 



1 Introduction 

The concept of noise in Physics stems from the theory of 
Brownian motion in which the force exerted on a macro- 
scopic impurity by the surrounding molecules is repre- 
sented by a random function T{t) 0. This model pro- 
vides a microscopic explanation of the law of diffusion 
and yields the celebrated fluctuation-dissipation relation 
0. Assuming that the Brownian particle of position x(t) 
and velocity v(t) is also submitted to a harmonic poten- 
tial U (x) = |u 2 i 2 , the dynamic equation associated with 
this system is the linear Langevin equation: 



dt 



v{t) + jv{t)+uj 2 x{t) =T{t), 



(1) 



where 7 is the effective damping coefficient. 

The interplay between noise and nonlinearity genera- 
tes many original phenomena that make the study of ran- 
dom dynamical systems a rich and fascinating held. Noise 
in a nonlinear system may induce non-equilibrium phase 
transitions or improve the performance of a device 
via stochastic resonance || . Random fluctuations may be 
rectified into a directed motion when a particle in a fluc- 
tuating environment is submitted to a "ratchet-like" po- 
tential f|. 

A simple nonlinear noisy dynamical system is obtained 
by considering Eq. ([!]) with an anharmonic potential U (x) . 
If the random force is approximated by a Gaussian white 



noise, the stationary probability distribution function (PDF) 
of the energy, valid when t 3> 7 -1 , is given by the cano- 
nical distribution. To the best of our knowledge, the closed 
form of the distribution function of the nonlinear oscilla- 
tor's energy is not known when the correlation time r of 
the random force is non-zero (see [ ?| and references therein 
for approximate results valid in the small r limit). 

In this work, we study the motion of an undamped non- 
linear oscillator submitted to an additive noise. For Gaus- 
sian white noise (Section 0), we apply the method pre- 
sented in a recent article fqf, where we studied the dyna- 
mical behavior of an undamped nonlinear oscillator with 
a fluctuating frequency represented as a parametric noise. 
We calculate the probability distribution function of the 
energy in the long-time limit and match it to the distribu- 
tion obtained in presence of damping. The average energy, 
root-mean-square amplitude and velocity of the oscillator 
grow algebraically with time, with anomalous diffusion ex- 
ponents different from those obtained for parametric noise. 
We show in Section |^ that diffusion is reduced at large 
times when the additive noise has a non-zero correlation 
time. Anomalous diffusion exponents are calculated in a 
self-consistent manner for an Ornstein-Uhlenbeck random 
force. In a somewhat counter-intuitive fashion, the white- 
noise behavior is observed for sufficiently short times. The 
crossover from white to colored noise can be embodied in 
a single scaling function as shown by dimensional analysis. 
All the results obtained for a nonlinear oscillator submit- 
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ted to additive or multiplicative, white or colored Gaus- 
sian noise are summarized in Table 1. 



2 The nonlinear oscillator with additive white 
noise 

We study a particle trapped in a confining potential 11 (x) 
and subject to additive noise. We neglect dissipative ef- 
fects unless stated otherwise. The dynamics of the nonlin- 
ear oscillator is given by 



dt 2 



x(t) 



dU{x) 
dx 



•£(*), 



(2) 



where x(t) represents the position of the oscillator at time 
t. In this Section, the random noise £(t) is a Gaussian 
white noise of zero mean- value and of amplitude T>: 



m = o, 



(3) 



For analytical as well as numerical calculations, we 
shall interpret stochastic differential equations such as (||) 
according to Stratonovich's convention. Indeed, Strato- 
novich calculus appears naturally in the limit of a vani- 
shingly small correlation time It must be emphasized 
that this convention (as opposed to Ito calculus) intro- 
duces correlations between the noise and the dynamical 
variables at the same time t ||. Moreover, the mechanical 
oscillator defined by Equations (@,^|) is not conservative: 
noise feeds energy into the system (see Eqs. (15) and (|35|)). 
In the absence of damping, the oscillator's amplitude in- 
creases with time. 

For the potential U(x) to be confining, we must have 
U — > +oo when |x| — > oo. We shall restrict our analysis 
to the case where IA (x) is a polynomial, and in order to 
respect the x — > — x symmetry, U{x) is taken to be even 
in x. Hence, when |x| — > oo, 



x 2n 

U with n > 2 , 

2n 



(4) 



where the coefficient of x 2n has been set to l/(2n) after 
a suitable rescaling. As the amplitude x of the oscillator 
grows at large times, only the asymptotic behavior of U(x) 
for |x| — * oo is relevant and thus Eq. (H) reduces to 



dt 2 



x(t) + x(t) 



2n— 1 



£(*)• 



(5) 



The method we shall use to study Eqs. (Q^) is akin 
to that presented in our previous work on non line ar oscil- 
lators with multiplicative noise || . In Section |2.l| , we use 
the integrability properties of the deterministic nonlinear 
oscillator to write exact stochastic differential equations 
in energy-angle vari able s. We then derive equipartition 
relations in Section 2.2. Averaging out the fast angular 
variable we calculate the energy's probability dis- 

tribution function in the long-time limit (Section 2.2 ). In 



Section 2.4, we obtain the associated anomalous diffusion 



exponents and constants and study their dependence on 
the stiffness of the confining potential at infinity. Our nu- 
merical simulations are based on a time discretization pre- 
sented in 1l2|l , and described in detail in 101 . 



2.1 Energy-angle coordinates 

We rewrite the second-order stochastic differential equa- 
tion (JsJ) as a first order system in energy (or action) and 
angle variables. Without noise, the deterministic version 
ofEq. (|): 

x + x 2 "- 1 = , (6) 

is integrable because of energy conservation, where the 
energy is defined by 



E = -x 2 + —x 2n . 
2 2n 

The action-angle variables (I, (f>) associated with Eq. 
are 



(7) 



7 = 4 



(2nB)Sn / 2n 

2E dx 

n 



[■1 

4(2" +1 n)^ / v 7 ! - u 2n du 

Jo 



oc E 2™ 



4> = \n 



Vl - U 2n 



(8) 
(9) 



where the angle variable <f> is defined modulo the oscilla- 
tion period AK n , with 



1 du 



Vl — a 



2)i 



(10) 



The solution of Eq. (g) can be parametrized as a func- 
tion of energy E and angle cj> 

x = E 1 ' 2n S n ^), (11) 
x = (2n)^E 1 ' 2 S' n (4>), (12) 

where the hyper-elliptic function S n is defined as 

' (2„)V2» du 



S n {(j>) = Y <-> (f) = y/n 



VT^j 2 
du 



V 2n 



(13) 



From this definition, we obtain the following relation be- 
tween S n and its derivative 5^ 

^ (l-«)'. (14) 
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Writing Eq. (||) in terms of energy and angle variables, 
we obtain the following stochastic dynamical system 



(a) 



1 S n {<t>) 



E 

4> = (2nE)^ 



r£(<). 



(2n)^ (2nEy. 
With the help of the auxiliary variable fl 

Eqs. @ and @ are written in the simpler form: 



C2 = nS' n {4>)m, 

fl \^ _ Snjfy 
(2n)*/ O 



(15) 
(16) 



(17) 

(18) 
(19) 



These equations have been derived without any hypothe- 
sis on the function £ and are rigorously equivalent to the 
original equation (^]). 

Although the method is identical, this set of equa- 
tions differs from that obtained for parametric noise (See 
Eqs. (27)- (28) in [§). As a consequence, additive and mul- 
tiplicative random noises lead to different long-time be- 
haviors of the observables of the system (See table 1). 



2.2 Equipartition relations 

Noting that S n and S' n are bounded functions of </>, we de- 
duce from Eqs. ( |is| ) and ( [T9| ) that fl is a diffusive variable 
that scales as t 1 *^ whereas is a fast variable that scales 
like i 3 / 2 " 1 / 2 ™ (these assertions will be justified rigorously 
in the next section). In the long time limit, it is therefore 
justified to average the dynamics over the variations of the 
fast variable 4>. Using Eqs. (|TTl)-(|r^), and assuming that 
4> is uniformly distributed over the interval [0, 4K n ] of a 
period, we find the following relations: 



(E) = 



1 



2n 

„ 2 «\ 



x 2 ) 



(20) 
(21) 



Figure | shows that Eqs. 

are indeed verified with 

excellent accuracy. 

Because the transformation of variables (|Il|)-([r2|) is 
independent of the nature of the noise £(i), the statis- 
tical identities (p0|)-( ]2l|) are the same as those obtained 
for parametric noise |p[. Such generalized equipartition 
relations are valid whenever the asymptotic probability 
distribution function is independent of the angle variable 
and is a function of the energy only. This is the case, in 
particular, of the canonical Boltzmann-Gibbs distribution: 
Eqs. (20)-(pT|) are valid at thermodynamic equilibrium, as 
can be shown using standard techniques of statistical me- 
chanics in the canonical ensemble. 




1000 



1000 



Fig. 1. Nonlinear oscillator with additive white noise: Eq. (^|) 
is integrated numerically for T> — 1, with a timestep St, and 
averaged over 10 4 realizations for n — 2 (2n — 1 = 3), 8t = 
5 1(T 4 ; n = 3 (2n - 1 = 5), St = 5 1(T 4 ; n = 4 (2n - 1 = 7), 
St = 2 1CT 4 . Fig. (a): the first equipartition ratio (E(tY) / (v(t) 2 ) 
is close to the theoretical value given in Eq. (j20|): 3/4 for 
n = 2; 2/3 for n = 3; 5/8 for n = 4. Fig. (b): the second 
equipartition ratio (v(t) 2 } / (x(t) 2n ) is close to 1 for n — 2,3,4. 



2.3 Calculation of the effective PDF 

The Fokker-Planck equation corresponding to the system 
<M 111) reads 



(n-l)/n 



d t P 




(22) 



d a (f(ct>)dnf(ct>)i 



where we have defined 

f(<j>)=nS , n (4>) and g{<t>) = S n (cf>) 



(23) 



This equation governs the dynamics of the probability dis- 
tribution function P*(i7, 4>). As we already noticed, the 
angle variable <f> varies rapidly as fl grows. We carry out 
the averaging of Eq. ( |22] ) under the hypothesis that the 
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Fig. 2. Nonlinear oscillator with additive, Gaussian white 
noise. Fig. (a): the ratio {E)/t is found close to T>/2 (V = 
1) for n = 2,3,4. Fig. (b): scaling behavior of (x(t) 2 }. 
Fig. (c): the ratio (x 2 )/^^" yields the following generalized 
diffusion constants D x n) : D {2) = 0.553(5), D ( x ] = 0.525(5), 
Z)£ 4 ' ) = 0.499(4), in excellent agreement with the predictions 
of Eqs. (|3rj)-(p8|). Numerical data is obtained from the same 
simulations as in Fig. [j]. 



probability density, Pt(f2, </>), becomes uniform in <f> when 
t — > oo. The averaged evolution equation of the reduced 
probability density Pt(J?) is then 



d t P = / 2 «>) 



V 



(24) 



where we have used / = (d^g)/n. The notation / 2 (</>) re- 
presents the mean value of the function f 2 over one period 
4K n of <p. This constant has been calculated explicitly in 
[p| and is given by 



/ 2 (</>) 



S'M) 2 = 



n 2 (2n)- 
n + 1 



(25) 



The Fokker-Planck equation (|24j ) is equivalent to an 
effective first-order stochastic differential equation for the 
slow variable fl. Defining an effective Gaussian white noise 
£(t) such that 



kl (*)£(*')> =f>5{t- t') with V = 



n 2 (2n)- 
n + 1 



-V, (26) 



n 



4|) the effective Langevin equation 
i(t). (27) 



we deduce from Eq. 
for Q 

V 1 

Yn ~n 

Thus Q can be reinterpreted as a Brownian variable in a 
logarithmic potential. 

The scale-invariant effective Fokker-Planck equation 
(E4j) can be solved explicitly. We find 



f2i 



-n 2 /(2Vt) 



r(^) (2i>ty- 



(28) 



where r is the Euler Gamma function. Using Eq. (|17|), we 
deduce from Eq. ( p^ ) the asymptotic probability distribu- 
tion function of the energy variable: 



Pt{E) 



1 



i 

E 



(2n) 



2Vt 



exp 



(2n) 



■E 



2Vt 



_ (29) 
Moreover, the P.D.F. (E9h can be matched to the cano- 
nical distribution. Suppose that some damping is present 
in the system. The dynamical equation (0) now reads 



—x(t)+~f-x(t)+x(t) 



2n-l 



dt 



dt 



(30) 



where 7 is the friction coefficient. The stationary distri- 
bution associated with the stochastic equation (|30| ) is: 



Pcan(-B) = 



1 



n+V 

In , 



1 /gTgy 

E\ V ) 



exp 



27SI 



(31) 



This PDF becomes the Boltzmann-Gibbs canonical mea- 
sure when supplemented with the fluctuation-dissipation 
relation pi. The crossover time t c from the asymptotic 
distribution function for the undamped oscillator to the 
Gibbs measure is found by matching Eq. ( |29| ) with Eq. (|T]) : 



tr 



2nj 



(32) 



Thus, there are three distinct time scales involved: a fast 
time scale over which the angle variable 4> gets uniformly 
distributed, a 'long' time scale over which the effective 
Fokker-Planck equation for i? ( p4] ) is relevant, and finally 
times longer than the crossover time t c beyond which the 
inevitable damping in the system takes over. 
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2.4 Time-asymptotic behavior of observables 



With the probability distribution function (|2S|), we can 
calculate the long time behavior, in the absence of damp- 
ing, of the expectation value of any observable that de- 
pends on the position and velocity of the particule. In 
particular, we obtain 



r 



( 2 ») 
r(4) 

nV 



2n 2 



1 



PI I , (33) 

(34) 
(35) 



From Eq. (p3[), we notice that the particle diffuses more 
slowly in a stiffer potential well, as one would intuitively 
expect. The equipartition relation between energy and ve- 
locity obtained in Section 2.2 is confirmed by Eqs. ( |34] , |35| ). 
For all potentials growing algebraically at infinity, we find 
that the mean value of the energy grows linearly with time, 
with a diffusion constant equal to V/2 (Eq. (p5|)). Such a 
remarkable and simple behavior was not obvious a priori. 

These predictions are confirmed by numerical simula- 
tions. The universal behavior of (E) is represented in Fi- 
gure ^.a. Using a numerical value of the noise amplitude 
T> = 1 in Eq. (B3), we obtain 



For n = 2, 
For n = 3, 
For n = 4, 



(x 2 ) = 0.552 f 1/2 



(x 2 ) = 0.526 f 1/3 , 
(x 2 ) = 0.500 f 1/4 . 



(36) 
(37) 
(38) 



These formulae are in excellent agreement with the nu- 
merical results displayed in Figures §.b and g.c. 

The validity of Eqs. ( p3p^ , |35| ) can be checked in the 
linear case, n = 1. The exact solution of the equation 



d 2 

— X(t)+L0 2 x(t)=at), 

is given by 

1 f* 

x(t) = — sin (uj(t — u)) du , 
w Jo 

x(t) = / cos(w(f — u)) £(u)du. 
Jo 

We obtain the following exact formulae (2j : 



rf 1 



V 

2J 2 



(i 2 > = ?* ( 1 



(E) 



2 
V 



sin 2wf 
2ut 
sin 2wf 
2uot 



-t . 



(39) 

(40) 
(41) 

(42) 
(43) 
(44) 



In the long time limit t — > oo, these expressions agree with 
Eqs. (|33| , p5| ) when one substitutes n = 1 (and w = 1). 



3 The nonlinear oscillator with additive 
colored noise 

In order to model an additive Gaussian noise with a non 
vanishing correlation time r, we use an Ornstein-Uhlenbeck 
process, denoted by rj(t). We thus replace Eq. (||) by the 
following equation 



dt 2 



x(t) + x(t) 



2n-l 



where the random noise 77(f) is obtained from 



dr?(f) 
df 



T T 



(45) 



(46) 



and £(f) is the Gaussian white noise defined in Eq. (||). 
The stationary statistical properties of 77 are given by 



(»?(*)) =0, 

l T 



(47) 



We rewrite Eq. (f 
as above [see Eqs. (G" 



(a) in (J7, <fi) coordinates and obtain 
f and @] 



r2 = n5;^)rj(t), 

12 



(48) 
(49) 



Again, the angle is a fast variable: we expect that the 
equipartition relations (^0|) and ( pl| ) are also valid for co- 
lored additive noise. This is indeed confirmed by numerical 
simulations (See F igure H ) . 

The system ( [^3|j4^j49| ) is equivalent to a Fokker-Planck 
equation for the P.D.F. Pj(/2, <f>, 77). However, averaging 
this equation over <fi does not lead to any conclusive re- 
sult. The qualitative reason for this failure is as follows: 
when the period of the angular variable becomes smaller 
than the coherence time r of the noise, the noise itself is 
averaged out. 

The goal of this section is to introduce alternative 
methods to determine the scaling behavior of nonlinear os- 
cillators with colored additive noise, first qualitatively in 
Section 3T thanks to dimen sion al analysis, then in a more 
rigorous manner in Section 3.2 usin g a self-consistent ar- 
gument. Finally, we show in Section 3.3 that the crossover 



between the behaviors typical of white noise and colored 
noise can be embodied in a single scaling function. 



3.1 Scaling analysis 



The equations ([15|) and ( |46| ) can be interpreted as a system 
of three coupled first-order stochastic differential equa- 
tions in presence of white noise. The corresponding Fokker- 
Planck equation for the P.D.F. P t (x,v,rj) is given by 



dP 
~dt 



-'/) 



dP 1 dnP V d 2 P 
\ 1 

dv t drj 2t 2 drf 



(50) 
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(a) 




1000 



1000 



Fig. 3. Nonlinear oscillator with additive, colored noise: 
Eqs. (p3|)-(|46|) are integrated numerically for T> = 1 and 
t = 1, with a timestep St = 10~ 6 and averaged over 
500 realizations for n = 2,3,4. Fig. (a): the first equiparti- 
tion ratio (E(tY}/(v(t) 2 } is close to the theoretical value ^1 
given in Eq. (EOJ). Fig. (b): the second equipartition ratio 
(v{t) 2 )/(x(t) 2n ) is close to 1. 



If we perform a scaling analysis in the spirit of |D| by com- 
paring terms in this equation two by two, we find that a 
consistent balance is found if v 2 ~ x 2n , a; 2 ™" 1 ~ r\ and 
rf ~ Vt/2 r 2 . Eliminating 77, we obtain the following scal- 
ing laws 



E 



Vt 
2^2 



Vt 
2^2 



Vt 



(51) 



A similar argument was used in |g] in the case of colored 
multiplicative noise, and yielded different exponents (See 
Table 1). 

Numerical simulations confirm the scalings predicted 
by Eqs. ( |5l"| ) for general n (See Figure ||). Besides, these 
scalings are consistent with exact results found in the li- 
near case n = 1 . In the next section we present an analytic 
derivation of the exponents appearing in Eq. (|5l| ) . 



3.2 A self-consistent calculation 

We make the a priori Ansatz that, in the long-time limit, 
Q grows algebraically with time as 



n ~ f 



(52) 



and we will determine the scaling exponent a from a self- 
consistent calculation. We deduce from Eq. ([l9]) that <f> ~ t v 
where 

n — 1 , . 

v= a+1. (53) 



Substituting this scaling of 4> in Eq. (flq), we obtain (here- 
after we shall leave aside all proportionality constants) 



(54) 



We must now determine the asymptotic statistical be- 
havior of this expression in the t — > 00 limit. In order to 
simplify the calculation, we replace the Ornstcin-Uhlcnbeck 
noise by a discrete dichotomous noise, where time is dis- 
cretized in steps of duration r (see [[l4| for a general dis- 
cussion). This approximation amounts to replacing expo- 
nentially decaying correlations, by finite time correlations 
and therefore leaves the diffusion exponents unchanged. 
Equation ( |54|) then reduces to 



/ dzS£(z* 

fc=0 



(55) 



kr 



where the variable takes the value ±^V/(2t) randomly 
during the time interval [kr , (k + l)r]; in other terms, 



(efe£;) = TT^kl 
At 



From Eqs. 



(n 2 




5|) and (|56j), we deduce that 

c(fe+l)T 

k=o \ Jkj 
Integrating by parts we obtain 



dzS'(z" 



(56) 



(57) 



(fe+l)T 



dz s^n - 

S n {k v T v 



S n ((k+l)<v) 

i/((k- 



-l)r) 

I /-(fc+l)T 



>(k T y 



dz 



5 n (z") 



(58) 



The integral term on the r.h.s. of (|58| ) is of order 0{k~ v ). 
Since S n is a bounded function, the first two terms are 
of order (D(k~" +1 ), and will dominate when k ^ 1 since 
v > 1. Wc thus obtain 



{Q 2 



1 



t/T 

k=i I kr] 



(59) 



Kirone Mallick, Philippe Marcq: Scaling behavior of a nonlinear oscillator with additive noise, white and colored 



7 



(a) 




1000 



(h) 




1000 



Fig. 4. Nonlinear oscillator with additive Ornstein-Uhlenbeck 
noise. The ratios (v 2 ) /t?^ (pig. (a)) and (x 2 )/t^^ 
(Fig. (b)) are approximately constant at large time for n = 
2, 3, 4. Numerical data is obtained from the same simulations 
as in Fig. |[ 

Assuming that the variable fl is not multifractal, we 
know from the scaling hypothesis, Eq. (|52|), that J? 2 ~ t 2a 



Comparing with Eq. (|59|), we deduce that 
2a = 3 - 2v . 



(60) 



Equations ( |53|) and ( |60| ) provide the required self-consistent 
condition and allow to calculate the exponents: 



2(2n- 1) 



and v = 



5n — 3 
2(2n- 1) 



(61) 



We can now deduce from Eq. ( |l7|) and from the equipar- 
tition relations, Eqs. ( |2C| ) and (|21|), the scaling behavior 
of the dynamical variables 



v 



t— , 

tat 2 "- 1 ) 
1 



(62) 

This result is in agreement with the scaling analysis of the 
previous section [see Equation (|5l|)]. 



3.3 Crossover from white to colored noise 

We shall now match the white and the colored noise be- 
haviors obtained in the long-time limit through a sca- 



ling function and will describe the crossover from the first 
regime to the second one. In the system ( f45| ) with Ornstein- 
Uhlenbeck noise, one can define two dimcnsionless vari- 
ables: t/r and T>r . Any dimensionless quantity can 
be written as a function of these two variables (or of any 
two independent combinations of these two variables). For 
example, the mean value of the energy can be written in 
full generality, as a function of time, as 



in 3n-l 

Vr~t, 



(63) 



where the prefactor Vt/2 ensures that the equation is di- 
mensionally correct and has been chosen to match the 
white noise limit. The scaling function <P allows to inter- 
polate between white and colored noise and exhibits the 
following asymptotic behavior: 

— The white noise limit corresponds to r = 0, therefore 
<Z>(0,0) = 1. 

— When r is finite (colored noise case) and t — * 00, wc 
deduce from Eq. (|5l| ) that, for u — > 00 and for v finite, 

1 — n 

<P(u,v) — > u 2 ™- 1 4>{v). The prefactor <ft(v) is indepen- 
dent of time and cannot be determined by dimensional 
analysis. 

The crossover between white noise and colored noise is 

2n 

obtained when u = Dr™- 1 1 ~ 1 , i.e., for a typical time of 
the order of 

_ 271 

T 1-1 

Across ^ ^ • (6^0 

At this crossover time the energy is of the order of E ~ 
~ t ™ _1 • The period of the deterministic oscilla- 
tor for such an energy is given by T ~ E~^~ ~ r (see 
Eq. (|l6|)). Thus, the crossover occurs when the period of 
the deterministic oscillator is of the order of the correla- 
tion time of the noise. This fact has the following intuitive 
explanation. When t -C Across, the period of the determi- 
nistic oscillator is very large compared to the correlation 
time of the noise: the noise is totally uncorrelated over 
one period and can be viewed as a sequence of uncorre- 
lated kicks and thus acts as a white noise. However, when 
t 3> icross, the period of the deterministic oscillator is small 
compared to the correlation time of the noise: the noise is 
highly correlated over a period and can be viewed as an 
almost constant quantity over this time scale. Many peri- 
ods must be added before the effect of the noise becomes 
perceptible: hence the diffusion slows down. 

We remark that in Eq. (|63| ) we have not used t/r as 
the dimensionless time variable which would have been the 
most natural choice. In fact, the white noise limit r — > 0, 
and the long time limit t — > 00 with finite r are indistin- 
guishable if one uses the variable t/r. For this reason, it 
may be incorrectly stated that 'in the long time limit, co- 
lored noise appears as white'. We have seen in our problem 
that exactly the opposite is true: the effect of a non-zero 
correlation time becomes relevant at long times. The limit 
t/r — * 00 is singular and its value depends on whether 

3n-l 

the second dimensionless variable v = T>r is equal 
to zero or is strictly positive. By choosing the variable 
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GAUSSIAN NOISE 



ADDITIVE 



MULTIPLICATIVE 



WHITE 



X ~ t 2 " 
1 

X ~ t 2 



p 



a; ~ t 2 <..-i) 



E ~ if"- 1 ) 



COLORED 



£2(2 


i-l) 


t 2 (2 


a-1) 




i-l) 


1 




-1) 


t 4(r 


-1) 


> *=>(' 


>-l) 



Table 1. Anomalous diffusion exponents for an undamped 
noisy nonlinear oscillator. 



u = T>T~^it as the dimensionless time variable in Eq. ( p3[ ) 
rather than t/r we avoid this difficulty because the white 
noise and the colored noise cases now appear as different 
limits. 



In Table 1, we summarize and compare the results that 
we have derived for additive noise and multiplicative noise, 
in this work and in , respectively. In the case of white 
noise, precise results are available: thanks to the averaging 
technique, the asymptotic probability distribution func- 
tions are known. For colored noise, only the anomalous 
diffusion exponents have so far been calculated. We no- 
tice that in all cases the exponent for the velocity is half 
the exponent for the energy and n times the exponent for 
the amplitude i.e., the ratios between corresponding expo- 
nents are independent of the problem considered. This is 
the consequence of the generalized equipartition relations, 
which are independent of the nature of the noise (in fact, 
equipartition also provides universal relations between the 
generalized diffusion constants appearing asprefactors in 
the scaling laws). Besides, it was argued in M that, in the 
long time limit, the multiplicative noise always dominates 
over the additive noise and indeed we observe that the 
diffusion exponents for the multiplicative noise are always 
larger than those for the additive noise. 

Although we have been able here to derive the ex- 
ponents for additive colored noise, the generalized diffu- 
sion constants, as well as the long time asymptotics of the 
probability distribution function were not found by our 
approach and their calculation remains an open problem. 
Besides, the precise calculation of the scaling function that 
describes the crossover between white and colored noise is 
certainly a challenging problem. 



4 Conclusion 

A Hamiltonian nonlinear oscillator submitted to random 
internal noise (additive noise) exhibits anomalous diffu- 
sion phenomena. We have calculated analytically the as- 
sociated scaling exponents in the case of white noise and 
have obtained, in the long time limit, an explicit expres- 
sion for the probability distribution function in phase spa- 
ce. The energy of such an oscillator grows linearly with 
time irrespective of the form of the confining potential. It 
is remarkable that even the rate of energy growth does not 
depend on the stiffness of the potential and is simply pro- 
portional to the amplitude of the noise. Our results also 
describe the behavior of an oscillator with small damping 
rate 7 for times less than the dissipation time 7 -1 . We 
have shown that the probability distribution function of 
the energy of the undamped oscillator can be matched 
with the canonical Boltzmann-Gibbs distribution when 

In the case of colored noise, the scaling exponents are 
modified. Their values have been determined by a self- 
consistent calculation and the predictions of dimensional 
analysis have been confirmed. The diffusion exponents are 
reduced because the coherence of the noise over a period 
of the system makes the energy transfer less efficient. The 
growth of energy is slower than linear and in the limit of 
a very stiff potential, the energy grows only as the square- 
root of time. 



We would like to thank Jean-Marc Luck for his constant in- 
terest in our work, as well as Michel Bauer, Peter Holdsworth 
and Kunihiko Kaneko for useful comments. 
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